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MODEL OF RECONNECTION OF WEAKLY STOCHASTIC MAGNETIC 

FIELD AND ITS IMPLICATIONS 

A. Lazarian, 1 and E. Vishniac 2 
RESUMEN 

Favor de proporcionar un resumen en espanol. If you are unable to translate your abstract 
into Spanish, the editors will do it for you. We discuss the model of magnetic field reconnection in 
the presence of turbulence introduced by us ten years ago. The model does not require any plasma effects 
to be involved in order to make the reconnection fast. In fact, it shows that the degree of magnetic field 
stochasticity controls the reconnection. The turbulence in the model is assumed to be subAlfvenic, with the 
magnetic field only slightly perturbed. This ensures that the reconnection happens in generic astrophysical 
environments and the model does not appeal to any unphysical concepts, similar to the turbulent magnetic 
diffusivity concept, which is employed in the kinematic magnetic dynamo. The interest to that model has 
recently increased due to successful numerical testings of the model predictions. In view of this, we discuss 
implications of the model, including the first-order Fermi acceleration of cosmic rays, that the model naturally 
entails, bursts of reconnection, that can be associated with Solar flares, as well as, removal of magnetic flux 
during star-formation. 

ABSTRACT 

We discuss the model of magnetic field reconnection in the presence of turbulence introduced by us ten years 
ago. The model docs not require any plasma effects to be involved in order to make the reconnection fast. In 
fact, it shows that the degree of magnetic field stochasticity controls the reconnection. The turbulence in the 
model is assumed to be subAlfvenic, with the magnetic field only slightly perturbed. This ensures that the 
reconnection happens in generic astrophysical environments and the model does not appeal to any unphysical 
concepts, similar to the turbulent magnetic diffusivity concept, which is employed in the kinematic magnetic 
dynamo. The interest to that model has recently increased due to successful numerical testings of the model 
predictions. In view of this, we discuss implications of the model, including the first-order Fermi acceleration of 
cosmic rays, that the model naturally entails, bursts of reconnection, that can be associated with Solar flares, 
as well as, removal of magnetic flux during star-formation. 
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1. PROBLEM OF ASTROPHYSICAL 
RECONNECTION 

Plasma conductivity is high for most astrophys- 
ical circumstances. This suggests that "flux freez- 
ing" , where magnetic field lines move with the local 
fluid elements, is a good approximation within as- 
trophysical magnctohydrodynamics (MHD). 

What happens when magnetic field lines inter- 
sect? This is the central question of the theory of 
magnetic reconnection. In fact, the whole dynam- 
ics of magnetized fluids and the back-reaction of the 
magnetic field depends on the answer. 
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2. THE SWEET-PARKER MODEL VERSUS 
PETSCHEK MODEL 

The literature on magnetic reconnection is rich 
and vast (see e.g., Priest & Forbes 2000 and ref- 
erences therein). We start by discussing a robust 
scheme proposed by Sweet and Parker (Parker 1957; 
Sweet 1958). In this scheme oppositely directed mag- 
netic fields are brought into contact over a region of 
L x size (see Fig. 1). The diffusion of magnetic field 
takes place over the vertical scale A which is related 
to the Ohmic diffusivity by rj w V r A, where V r is the 
velocity at which magnetic field lines can get into 
contact with each other and reconnect. Given some 
fixed rj one may naively hope to obtain fast reconnec- 
tion by decreasing A. However, this is not possible. 
An additional constraint posed by mass conservation 
must be satisfied. The plasma initially entrained on 
the magnetic field lines must be removed from the 
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reconnection zone. In the Sweet-Parker scheme this 
means a bulk outflow through a layer with a thick- 
ness of A. In other words, the entrained mass must 
be ejected, i.e., pV r L x = p'VaA, where it is assumed 
that the outflow occurs at the Alfven velocity. Ig- 
noring the effects of compressibility, then p = p' and 
the resulting reconnection velocity allowed by Ohmic 

— 1/2 

diffusivity and the mass constraint is V r ~ VaTZ-l , 
where Tl^ 1 ^ 2 = (ti/VaLx) 1 ^ 2 is the Lundquist num- 
ber. This is a very large number in astrophysical 
contexts (as large as 10 20 for the Galaxy) so that 
the Sweet-Parker reconnection rate is negligible. 

One of the ways to increase the reconnection rate 
is to change the global geometry thereby reducing 
L x . An example of the latter is the suggestion by 
Petschek (1964) that reconnecting magnetic fields 
would tend to form structures whose typical size in 
all directions is determined by the resistivity ('X- 
point' reconnection). This results in a reconnection 
speed of order Va/ ItiTZl- However, attempts to pro- 
duce such structures in simulations of reconnection 
have been disappointing (Biskamp 1984, 1986). In 
numerical simulations the X-point region tends to 
collapse towards the Sweet-Parker geometry as the 
Lundquist number becomes large (Biskamp 1996). A 
general review of astrophysical magnetic reconnec- 
tion theory can be found in Priest & Forbes (2003). 

Pctschek-type reconnection can survive in the 
presence of localized increase of effective resistivity, 
which can happen due to various plasma effects. For 
instance, recent years have been marked by a sub- 
stantial progress in simulations of collisionless recon- 
nection (see Shay & Drake 1998, Bhattacharjee et 
al. 2003, Drake et al. 2006a and references therein). 
This work indicates that under some circumstances 
a kind of standing whistler mode can stabilize an 
X-point reconnection region. However, these studies 
have not demonstrated the possibility of fast recon- 
nection for generic field geometries. They assume 
that there are no bulk forces acting to produce a 
large scale current sheet, and that the magnetized 
regions are convex, which minimizes the energy re- 
quired to spread the field lines. In addition, the 
requirements on the media being collisionless are 
rather strict, for instance, estimates in the litera- 
ture (see Uzdensky 2006, Yamada et al. 2006) sug- 
gest that the L x should not exceed approximately 50 
electron mean free paths. This makes the model not 
applicable to many astrophysical environments, e.g. 
to the interstellar medium. 

In any case, while the researchers argue whether 
Hall MHD or fully kinetic description (Daughton 
2006) is necessary, one statement is definitely true. 



If magnetic reconnection is only fast in collisionless 
environments, most of the MHD simulations, e.g. of 
interstellar medium, accretion disks, stars, where the 
environment is collisional, are in error. We shall ar- 
gue below that this radical conclusion may not be 
true and the reconnection is also fast in most astro- 
physical collisional environments. 

3. EFFECTS OF MAGNETIC STOCHASTICITY 
ON RECONNECTION 
3.1. Earlier Attempts 

The notion that magnetic field stochasticity 
might affect current sheet structures is not unprece- 
dented. In earlier work Speiser (1970) showed that 
in collisionless plasmas the electron collision time 
should be replaced with the time a typical elec- 
tron is retained in the current sheet. Also Jacob- 
son & Moses (1984) proposed that current diffusivity 
should be modified to include diffusion of electrons 
across the mean field due to small scale stochasticity. 
These effects will usually be small compared to effect 
of a broad outflow zone containing both plasma and 
ejected shared magnetic flux. Moreover, while both 
of these effects will affect reconnection rates, they are 
not sufficient to produce reconnection speeds com- 
parable to the Alfven speed in most astrophysical 
environments. 

" Hyper- resistivity" (Strauss 1985, Bhattacharjee 
& Hameiri 1986, Hameiri & Bhrattacharjee 1987, Di- 
amond & Malkov 2003) is a more subtle attempt to 
derive fast reconnection from turbulence within the 
context of mean-field resistive MHD. The form of the 
parallel electric field can be derived from magnetic 
hclicity conservation. Integrating by parts one ob- 
tains a term which looks like an effective resistivity 
proportional to the magnetic hclicity current. There 
are several assumptions implicit in this derivation. 
The most important objection to this approach is 
that by adopting a mean-field approximation one is 
already assuming some sort of small-scale smearing 
effect, equivalent to fast reconnection. Furthermore, 
the integration by parts involves assuming a large 
scale magnetic helicity flux through the boundaries 
of precisely the form required to drive fast reconnec- 
tion. Straus (1988) partially circumvented the first 
problem by examining the effect of tearing mode in- 
stabilities within current sheets. However, the result- 
ing reconnection speed enhancement is roughly what 
one would expect based simply on the broadening of 
the current sheets due to internal mixing. This ef- 
fect does not allow us to evade the constraints on 
the global plasma flow that lead to slow reconnec- 
tion speeds, a point which has been demonstrated 
numerically (see Matthaeus & Lamkin 1985). 
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Fig. 1. Upper plot: Sweet-Parker model of reconnection. 
The outflow is limited by a thin slot A, which is de- 
termined by Ohmic diffusivity. The other scale is an 
astrophysical scale L ^> A. Middle plot: Turbulent re- 
connection model that accounts for the stochasticity of 
magnetic field lines. The outflow is limited by the diffu- 
sion of magnetic field lines, which depends on field line 
stochasticity. Low plot: An individual small scale re- 
connection region. The reconnection over small patches 
of magnetic field determines the local reconnection rate. 
The global reconnection rate is substantially larger as 
many independent patches come together. 

3.2. Model in Lazarian & Vishniac 99 

The model of reconnection proposed in Lazarian 
& Vishniac (1999, henceforth LV99) is a natural gen- 
eralization of the Sweet-Parker model (see Fig. [T|). 
The problem of the Sweet-Parker model is that the 
reconnection is negligibly slow for any realistic as- 
trophysical conditions. However, astrophysical mag- 
netic fields are generically turbulent. 

LV99 consider the case in which there exists a 
large scale, well-ordered magnetic field, of the kind 
that is normally used as a starting point for discus- 
sions of reconnection. In addition, we expect that 
the field has some small scale 'wandering' of the field 
lines. On any given scale the typical angle by which 
field lines differ from their neighbors is <j> -C 1, and 
this angle persists for a distance along the field lines 
Aii with a correlation distance Aj_ across field lines 
(see Fig.HJ. 

To do quantiative estimates one has to adopt 
the model of MHD turbulence. Such a model for 



realistic compressible fluids may be constructed on 
the basis of the Goldrcich-Sridhar (1995, henceforth 
GS95) model of incompressible turbulence (see Cho 
& Vishniac 2000, Maron & Goldrcich 2001, Lithwick 
k Goldrech 2001, Cho, Lazarian & Vishniac 2002, 
Cho & Lazarian 2002, 2003). The most important 
in terms of magnetic reconnection are the properties 
of Alfvenic modes that determine field wandering as 
shown in LV99. 

The modification of the global constraint induced 
by mass conservation in the presence of a stochas- 
tic magnetic field component is self-evident. Instead 
of being squeezed from a layer whose width is de- 
termined by Ohmic diffusion, the plasma may dif- 
fuse through a much broader layer, L y ~ (y 2 ) 1 ^ 2 
(see Fig. [TJ , determined by the diffusion of magnetic 
field lines. This suggests an upper limit on the re- 
connection speed of ~ VaUd 2 ) 1 ^ 2 /L x ). This will be 
the actual speed of reconnection if the progress of 
reconnection in the current sheet does not impose a 
smaller limit. The value of (y 2 ) 1 ^ 2 can be determined 
once a particular model of turbulence is adopted, but 
it is obvious from the very beginning that this value 
is determined by field wandering rather than Ohmic 
diffusion as in the Sweet-Parker case. 

What about limits on the speed of reconnection 
that arise from considering the structure of the cur- 
rent sheet? In the presence of a stochastic field 
component, magnetic reconnection dissipates field 
lines not over their entire length ~ L x but only 
over a scale An <C L x (see Fig. [I}, which is the 
scale over which magnetic field line deviates from 
its original direction by the thickness of the Ohmic 
diffusion layer AJ 1 ps f)/V reCj i oca i. If the angle <fi 
of field deviation does not depend on the scale, the 
local reconnection velocity would be ~ Va4> an d 
would not depend on resistivity. In LV99 it is taken 
into account that <fi does depend on scale. There- 
fore the local reconnection rate V rec ^ oca i is given 
by the usual Sweet-Parker formula but with An in- 
stead Of L x , i.e. Vrecjocal W V A ( V A A|| / rj) ~ 1/2 . It 

is obvious from Fig. [1] that ~ L x /Xu magnetic field 
lines will undergo reconnection simultaneously (com- 
pared to a one by one line reconnection process for 
the Sweet-Parker scheme). Therefore the overall 
reconnection rate may be as large as V reCtg i b a i ~ 
V A {L X / \\){V A \\/iiy 1/2 ■ Whether or not this limit 
is important depends on the value of A|| . 

The relevant values of A|| and {y 2 } 1 ^ 2 depend on 
the magnetic field statistics. This calculation was 
performed in LV99 using the GS95 model of MHD 
turbulence providing the upper limit on the recon- 
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nection speed: 



V, 



r,up 



Va min 



L, 



where I and vi are the energy injection scale and 
turbulent velocity at this scale respectively. In LV99 
other processes that can impede reconnection were 
found to be less restrictive. For instance, the tan- 
gle of reconnection field lines crossing the current 
sheet will need to reconnect repeatedly before indi- 
vidual flux elements can leave the current sheet be- 
hind. The rate at which this occurs can be estimated 
by assuming that it constitutes the real bottleneck 
in reconnection events, and then analyzing each flux 
element reconnection as part of a self-similar system 
of such events. This turns out not to impede the re- 
connection. As a result, LV99 concludes that Jl} is 
not only an upper limit, but is the best estimate of 
the speed of reconnection. The thick plasma outflows 
observed during the 2003 November 4 Coronal Mass 
Ejection reported in Ciaravella & Raymond (2008) 
are also consistent with the LV99 model. 

3.3. Generalizing the model for partially ionized gas 

A partially ionized plasma fills a substantial vol- 
ume within our galaxy and the earlier stages of star 
formation take place in a largely neutral medium. 
This motivates our study of the effect of neutrals 
on reconnection. The role of ion-neutral collisions 
is not trivial. On one hand, they may truncate the 
turbulent cascade, reducing the small scale stochas- 
ticity and decreasing the reconnection speed. On the 
other hand, the ability of neutrals to diffuse perpen- 
dicular to magnetic field lines allows for a broader 
particle outflow and enhances reconnection rates. 

Reconnection in partially ionized gases before 
the introduction of the LV99 model looked hope- 
lessly slow. For instance, in Vishniac & Lazarian 
(1999, henceforth VL99) we studied the diffusion of 
neutrals away from the reconnection zone assum- 
ing anti-parallel magnetic field lines The ambipo- 
lar reconnection rates obtained in VL99, although 
large compared with the Sweet-Parker model, are 
insufficient either for fast dynamo models or for 
the ejection of magnetic flux prior to star forma- 
tion. In fact, the increase in the reconnection speed 
stemmed entirely from the compression of ions in the 
current sheet, with the consequent enhancement of 
both recombination and ohmic dissipation. This ef- 
fect is small unless the reconnecting magnetic field 
lines are almost exactly anti-parallel (VL99, see also 
Heitsch & Zweibel 2003ab). Any dynamically sig- 



nificant shared field component will prevent notice- 
able plasma compression in the current sheet, and 
lead to speeds practically indistinguishable from the 
standard Sweet-Parker result. Since generic recon- 
nection regions will have a shared field component 
of the same order as the reversing component, the 
implication is that reconnection and ambipolar diffu- 
sion do not change reconnection speed significantly. 

Lazarian, Vishniac & Cho (2004, henceforth 
LVC04) presented a model of turbulence in a par- 
tially ionized gas. This model agrees well with nu- 
merical simulations available as the limiting case 
which can be characterized by one fluid with a high 
Prandtl number (see Cho et al. 2002, 2003). Using 
this model LVC04 described field wandering, which 
is the core of the LV99 model of reconnection. They 
showed that the magnetic reconnection proceeds fast 
both in diffuse interstellar and molecular cloud par- 
tially ionized gas. 

4. FIRST ORDER FERMI ACCELERATION 
INDUCED BY STOCHASTIC 
RECONNECTION 

In Sweet-Parker model the reconnection can ac- 
celerate charged particles, e.g. due to the elec- 
tric field in the reconnection region. However, the 
speed of Sweet-Parker reconnection is negligible for 
most astrophysical environments, thus the transfer 
of energy from the magnetic field to particles is ab- 
solutely negligible, if the reconnection follows the 
Sweet-Parker predictions. 

It is interesting to notice that the first-order 
Fermi acceleration process is intrinsic to the LV99 
model of reconnection. To understand it consider 
a particle entrained on a reconnected magnetic field 
line (see FigfT]). This particle may bounce back and 
forth between magnetic mirrors formed by oppositely 
directed magnetic fluxes moving towards each other 
with the velocity Vr. Each of such bouncing will 
increase the energy of a particle in a way consistent 
with the requirements of the first-order Fermi pro- 
cess (de Gouveia Dal Pino & Lazarian 2001, 2003, 
2005, Lazarian 2005). 

Another way of understanding the acceleration of 
energetic particles in the reconnection process above 
is to take into account that the length of magnetic 
field lines is decreasing during reconnection. As a 
result, the physical volume of the energetic particles 
entrained on the field lines is shrinking. Thus, due to 
Louivillc theorem, their momentum should increase 
to preserve the constancy of the phase volume. 

An interesting property of this mechanism that 
potentially can be used to test the acceleration ob- 
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Fig. 2. Cosmic rays spiral about a reconnected mag- 
netic field line and bounce back at points A and B. The 
reconnected regions move towards each other with the 
reconnection velocity Vr. The advection of cosmic rays 
entrained on magnetic field lines happens at the outflow 
velocity, which is in most cases of the order of Va ■ Bounc- 
ing at points A and B happens because either of stream- 
ing instability or turbulence in the reconnection region 
(from Lazarian 2005). 

servationally is that the resulting spectrum of accel- 
erated particles is different from that arising from a 
shock. Gouveia Dal Pino & Lazarian (2003, 2005) 
used this mechanism of particle acceleration to ex- 
plain the synchrotron power-law spectrum arising 
from the flares of the microquasar GRS 1915+105. 
Note, that the mechanism acts in the Sweet-Parker 
scheme as well as in the scheme of turbulent recon- 
nection. However, in the former the rates of recon- 
nection and therefore the efficiency of acceleration 
are marginal in most cases. 

The mechanism is similar to the acceleration 
mechanism that was proposed later by Drake et al. 
(2006b). Drake et al. (2006b) considered the acceler- 
ation of electrons and, similarly, to the Matthaeus, 
Ambrosiano & Goldstein (1984), assumed that the 
acceleration happens within 2D contracting loops. 
While for LV99 model of reconnection the generic 
configuration of magnetic field arc contracting spirals 
this does not induce a radical difference between the 
processes. The difference in the spectrum of acceler- 
ated particles obtained by Drake et al. (2006b) and 
that obtained by us stems from the fact that in the 
former paper the strong backreaction of the accel- 
erated electrons was assumed. However, depending 
on the environment, this backreaction may vary, re- 
sulting in variations of the spectral slopes. Further 
studies of the acceleration process (see Fig. [2]) are 



necessary. 

5. FLARES INDUCED BY STOCHASTIC 
RECONNECTION 

If turbulence can drive reconnection, which in 
turn transforms magnetic energy into kinetic energy, 
then it seems appropriate to wonder if the process 
can be self-sustaining. That is, given a very small 
level of ambient turbulence, how likely is it that re- 
connection will speed up as it progresses, without 
any further input from the surrounding medium? A 
detailed examination of this is beyond the scope of 
current calculations, but we can clarify this with a 
simple physical model. 

Let's consider a reconnection region of length L 
and thickness A. The thickness is determined by the 
diffusion of field lines, which is in turn determined by 
the strength of the turbulence in the volume. Recon- 
nection will allow the magnetic field to relax, creat- 
ing a bulk flow. However, since stochastic reconnec- 
tion is expected to proceed unevenly, with large vari- 
ations in the current sheet, we can expect that some 
unknown fraction of this energy will be deposited 
inhomogeneously, generating waves and adding en- 
ergy to the local turbulent cascade. We take the 
plasma density to be approximately uniform so that 
the Alfven speed and the magnetic field strength are 
interchangeable. The nonlinear dissipation rate for 
waves is 



nonlinear 



2 „,2 



(2) 



where the first rate is the self-interaction rate for 
the waves and the second is the dissipation rate by 
the ambient turbulence (see Beresnyak & Lazarian 
2008ab) The important point here is that k± for the 
waves falls somewhere in the inertial range of the 
strong turbulence. Eddies at that wavenumber will 
disrupt the waves in one eddy turnover time, which 
is necessarily less than L/Va- The bulk of the wave 
energy will go into the tubulent cascade before es- 
caping from the reconnection zone. (This zone will 
radiate waves, for the same reason that turbulence 
in general radiates waves, but it will not significantly 
impact that energy budget of the reconnection re- 
gion.) 

We can therefore simplify our model for the en- 
ergy budget in the reconnection zone by assum- 
ing that some fraction e of the energy liberated by 
stochastic reconnection is fed into the local turbu- 
lent cascade. The evolution of the turbulent energy 
density per area is 

V A 



4- 

dt 



turb ) 



eVWr 



A v rec 



J turb 
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where the loss term covers both the local dissipa- 
tion of turbulent energy, and its advection out of 
the reconnection zone. Since V rec ~ v tur b and 
A ~ L(vturb/VA), we can rewrite this by defining 
Ma = vturb /Va and r = L/Va so that 



d_ 



eM A - M z 



If e is a constant then 



v turb 



V A e^ 



1 



1 



Ml 
e 



3 -2t/3 



1/2 



(4) 



(5) 



This implies that the reconnection rate rises to 
e 1 / 21 ^ is a time comparable to the ejection time 
from the reconnection region (~ L/Va)- Given that 
reconnection events in the solar corona seem to be 
episodic, with longer periods of quiescence, this im- 
plies that either e is very small, for example - de- 
pendent on the ratio of the thickness of the current 
sheet to A, or is a steep function of Ma- If it scales 
as M.a to some power greater than two then initial 
conditions dominate the early time evolution. 

An alternative route by which stochastic recon- 
nection might be self-sustaining would be in the con- 
text of a series of topological knots in the magnetic 
field, each of which is undergoing reconnection. Now 
the problem is sensitive to geometry. Let's assume 
that as each knot undergoes reconnection it releases 
a characteristic energy into a volume which has the 
same linear dimension as the distance to the next 
knot. The density of the energy input into this vol- 
ume is roughly eV^vturb/L, where e is the efficiency 
with which the magnetic energy is transformed into 
turbulent energy. We have 



eV 



2 ^turb 



v' 3 
Lk ' 



(6) 



where is the distance between knots and v' is the 
turbulent velocity created by the reconnection of the 
first knot. This process will proceed explosively if 
v' > v^rb or 

V%L k e > v 2 turb L. (7) 

This condition is almost trivial to fulfill. The bulk 
motions created by reconnection will unavoidably 
generate significant turbulence as they interact with 
their surrounding, so e should be of order unity. 
Moreover the length of any current sheet should be 
at most comparable to the distance to the nearest 
distinct magnetic knot. The implication is that each 
magnetic reconnection event will set off its neigh- 
bors, boosting their reconnection rates from v tur b, 
set by the environment, to c^-^Va^L^/ L) 1 / 2 (as long 



as this is less than Va)- The process will take a 
time comparable to L/v tur b to begin, but once ini- 
tiated will propagate through the medium with a 
speed comparable to speed of reconnection in the in- 
dividual knots. In a more realistic situation, the net 
effect will be a kind of modified sandpile model for 
magnetic reconnection in the solar corona and chro- 
mosphere. As the density of knots increases, and 
the energy available through magnetic reconnection 
increases, the chance of a successfully propagating 
reconnection front will increase. 



6. STOCHASTIC RECONNECTION AND 
REMOVAL OF MAGNETIC FLUX FROM 
MOLECULAR CLOUDS 

As we mentioned above, Sweet-Parker reconnec- 
tion is too slo\41, while collisionless Petschek recon- 
nection is not applicable to molecular clouds. At 
the same time, Shu et al. (2006) showed that mag- 
netic field is being removed from the star-forming 
cores faster than it is allowed by the standard am- 
bipolar diffusion scenario (see Tassis & Mouschovias 
2005). Shu et al. (2007) proposed a mechanism that 
required efficient reconnection of magnetic loops re- 
ferring to the "hyperresistivity" concept. As we dis- 
cussed earlier, this concept is not self-consistent and 
problematic at its core. Could the reconnection be 
done by the mechanism we discuss in the paper? 

As we discussed earlier, LVC04 considered mag- 
netic reconnection in partially ionized gas and ob- 
tained fast reconnection rates that can be obtained 
there. Thus, it is suggestive that magnetic field can 
be removed this way from molecular clouds. Inciden- 
tally, this also means that the model of reconnection 
should be considered for the transport of the angular 
momentum in protostellar disks (see Lazarian 2005). 

7. DISCUSSION AND SUMMARY 

The advantage of the model in Lazarian & Vis- 
haniac (1999) is that, first of all, the reconnection 
is robust and is fast in any type of fluid, provided 
that the fluid is turbulent enough. The latter re- 
quirement is natural for most of astrophysical fluids 
(see Armstrong et al. 1995). At the same time, the 
model predicts that if the fluids are not turbulent 
initially, they should be prone to bursts of reconnec- 
tion, which may provide an appealing explanation of 
Solar Flares. 



3 A modification of the Sweet-Parker model to include am- 
bipolar diffusion cannot generically induce much faster recon- 
nection either (see Vishniac &; Lazarian 1999). 
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7.1. Self- consistency of the model 

The high speed of reconnection given by equation 
(fT]) naturally leads to a question of self-consistency. 
Is it reasonable to take the turbulent cascade sug- 
gested in GS95 when field lines in adjacent eddies 
are capable of reconnecting? It turns out that in 
this context, our estimate for V reCl giobal is just fast 
enough to be interesting. We note that when con- 
sidering the intersection of nearly parallel field lines 
in adjacent eddies the acceleration of plasma from 
the reconnection layer due to the pressure gradient 
is not k\\V\, but rather (fcjj/fc^ )V^, since only the 
energy of the component of the magnetic field which 
is not shared is available to drive the outflow. On 
the other hand, the characteristic length contraction 
of a given field line due to reconnection between ad- 
jacent eddies is only k\\/k\. This gives an effective 
ejection rate of fcn Va- Since the width of the dif- 
fusion layer over a length k7 l is just fc7 , we can 



replace equation (fT]) with V rec , g i bai 
associated reconnection rate is just 



reconnect 



V A k h 



The 



(8) 



which in GS95 is just the nonlinear cascade rate on 
the scale k7 . However, this result is general and 
does not involve assuming that GS95 is correct. As 
we discuss below, most of the energy liberated in re- 
connection goes into motions on length scales com- 
parable to the dimensions of the reconnecting eddies, 
so this energy release will not short circuit the en- 
ergy cascade described in GS95. On the other hand, 
we can invert this argument to sec that reconnec- 
tion can play an important role in preventing the 
buildup of unresolved knots in the magnetic field. 
Such structures could play a major role in inhibiting 
the cascade of energy to smaller scales, flattening 
the energy spectrum relative to the predictions of 
GS95. Our conclusion is that such structures will 
disappear as fast as they appear, supporting the no- 
tion that they play a limited role in the dynamics of 
MHD turbulence. 

Finally, we note that if the magnetic field struc- 
ture is driven by turbulence in another location, as 
when the footpoints of magnetic arcades are stirred 
by turbulent motions, then we can evaluate its effects 
in terms of the amplitude of field stochasticity and 
the scaling of structure anisotropy with scale. The 
actual turbulence may be balanced or imbalanced, 
have or not have "polarization intermittency" (see 
Beresnyak & Lazarian 2006, 2008) the robust na- 
ture of the reconnection implies that the reconnec- 
tion will be sensitive to the amplitude of the induced 



field stochasticity, but not the details of the turbu- 
lent mixing process. 

7.2. Turbulent diffusivity and dynamo 

To enable sustainable dynamo action and, for ex- 
ample, generate a galactic magnetic field, it is neces- 
sary to reconnect and rearrange magnetic flux on 
a scale similar to a galactic disc thickness within 
roughly a galactic turnover time (~ 10 8 years). This 
implies that reconnection must occur at a substan- 
tial fraction of the Alfven velocity. The preceding 
arguments indicate that such reconnection velocities 
should be attainable if we allow for a realistic mag- 
netic field structure, one that includes both random 
and regular fields. This does solve one part of the 
problem of dynamo. The other part is related to 
magnetic helicity conservation (see Vishniac et al. 
2003). 

Interestingly enough, high turbulent diffusivity 
is also required for advecting heat in astrophysi- 
cal plasmas, e.g. in clusters of galaxies (see Lazar- 
ian 2006). Turbulent reconnection enables eddies to 
transfer heat in the way similar to the advection of 
heat by turbulence in unmagnetized fluid (see Cho 
ct al. 2003). 

7.3. Dissipation of energy 

The usual assumption for energy dissipation in 
reconnection is that some large fraction of the energy 
given up by the magnetic field, in this case ~ pV\L\, 
goes into heating the electrons. This is not the case 
here. Only a fraction, ~ l/(k\\L x ) of any flux cl- 
ement is annihilated by Ohmic heating within the 
reconnection zone. Over the entire course of the re- 
connection event the efficiency for electron heating 
is no greater than 



n/A 



V, 



rec.local 



l 



Vrec : global Vrec, global 



k\\L x 



(9) 



or, with GS95 scaling substituted (see LV99) 



VaI_ 

T) 



-2/5 



8/5 



4/5 



(10) 



The electron heating within the current sheet will 
not be uniform, due to the presence of turbulence, 
the intermittent presence of reconnected flux, and 
any collective effects we have neglected here. To the 
extent that these are important they will also lower 
the electron heating efficiency by broadening the re- 
connection layer. 
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7.4. Summary 
The results above can be summarized as follows 

• Recent successful numerical testing of the 
model in Lazarian & Vishniac (1999), presented in 
a companion paper by Kowal ct al. (2009, this vol- 
ume) increase the appeal of the model and stimulate 
studies of its consequencies. 

• The aforementioned model of stochastic field 
reconnection provides justification for many of as- 
trophysical simulations. If, on the contrary, the only 
reconnection model that works is the collisionless re- 
connection, this means that most of the numerical 
simulations, for instance, of interstellar medium are 
in error. Indeed, in the collisional environments the 
reconnection speed would be negligible, which can- 
not be achieved with the numerical simulations. 

• The model of magnetic field reconnection in the 
presence of weak magnetic field stochasticity is a nat- 
ural generalization of the Sweet-Parker model of re- 
connection. Its many consequencies include bursts of 
reconnection, efficient first order Fermi acceleration 
of energetic particles, efficient diffusion of magnetic 
field in turbulent plasmas etc. 
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